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ABSTRACT 

We investigate the impact of different initial conditions for the initial density profile 
and the initial turbulence on the formation process of protostellar clusters. We study 
the collapse of dense molecular cloud cores with three-dimensional adaptive mesh 
refinement simulations. We focus our discussion on the distribution of the gas among 
the protostellar objects in the turbulent dynamical cores. Despite the large variations 
in the initial configurations and the resulting gas and cluster morphology we find 
that all stellar clusters follow a very similar gas accretion behaviour. Once secondary 
protostars begin to form, the central region of a cluster is efficiently shielded from 
further accretion. Hence, objects located close to the centre are starved of material, 
as indicated by a strong decrease of the central accretion rate. This Fragmentation 
Induced Starvation occurs not only in rotationally supported discs and filaments, but 
also in more spherically symmetric clusters with complex chaotic motions. 
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1 INTRODUCTION 

The current paradigm of star formation suggests that most 



stars form in groups rather than in isolation (Lada et al 



|2003 ). However, the fraction of stars that form in a bound 



clustered environment is still a matter of debate and varies 



between different models (Bressert et al. 2010). Whereas 



massive stars only seem to form in a dense environment, 
low-mass stars form in all observed star-forming regions. Nu- 
merical simulations suggest the formation of massive stars 
simultaneously with the formation of a cluster fSmith et al 



2009[ [Peters et aL]|2010a|c|b| . The spatial distribution of 



stars in the forming cluster shows a dependence on the stef 
lar mass: more massive stars are located closer to the cen- 
tre of the group, low-mass stars tend to populate the outer 
regions. Whether this mass segregation is primordial and 
therefore determined by the formation scenario or due to 
a dynamical relaxation process is still a matter of debate 
and might differ significantly among different clusters (see, 
e.g., the reviews by [Mac Low fc Klcss cn] [20041 |McKee fc[ 



Scaio|[T986| [ScaIo|[T998l [Kroupa[[200l| [Chabrier[[2003l [BaF| 
tian et al.[[2010 |. Yet it is unclear, what is the influence of 



various physical processes and initial conditions on the ini- 
tial mass distribution in star forming regions. Different phys- 
ical processes like radiative feedback and magnetic fields are 
likely to have an impact on the collapse of gas clouds and the 
subsequent formation of a stellar cluster. Radiative feedback 
from protostars can have two opposite effects. On the one 
hand, it heats the surrounding gas and thus contributes to 
stabilising the gas against collapse by increasing the Jeans 
mass. In case of massive protostars, the strong contribution 
in UV ionises the gas, forming Hll regions. Hydrodynamic 
simulations of forming clusters show that radiative feedback 
tends to reduce the degree of fragmentation, but does not 



suppress it entirely (Krumholz et al. 2007 Bate 2009 Pe- 
ters et al.|2010a|c|b |. On the other hand, radiation pressure 



may indirectly enhance the formation of dense cores and 
the subsequent condensation into protostars. There are two 
different scenarios for this triggered star formation process. 
The first, "collect and collapse" (Elmegreen & Lada 1977 



in observations as well as in numerical simulations is a unf 
versal distribution of masses (initial mass function, IMF 
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gion pushes ambient gas into a shell, followed by fragmenta- 
tion and collapse of the swept-up material. The other one is 
called "radiation driven implosion" ( Bertoldi|1989 Bertoldi 



fc McKee|1990 1 and depicts the UV-driven compression of a 
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cloud core that is embedded into an Hll region. Apart from 
radiation, magnetic fields have a noticeable impact on the 
evolution of a collapsing core. Magnetic pressure and tension 
forces counter the gravitational collapse and diminish frag- 
mentation without completely inhibiting it ([Z icglcr 2 005[ 
Banerjee fc Pudritz|p006l [Price fc B ate''2007l Hennebelle 



fc Teyssier||2008| |Hennebelle fc Frm nang 2008 



Hennebelle 

fc Ciardi "2009': Commcrgon et al. '2010 Biirzlc et al. "2011[ 



Seifried et al..2011j). Likewise, different initial conditions of 



the gas cloud strongly influence the star formation process 



( Girichidis et al.|2011 hereafter Paper I). 

However, a clear universal formation picture of stellar 
systems in dense environments is still missing. Several for- 
mation and gas accretion processes have been proposed with 
different distribution scenarios of the gas onto the protostel- 
lar objects and different predictions concerning the more 
massive stars and the resulting IMF. The apparent similar- 
ity between the stellar mass function and the mass function 



of bound cores ( Testi & Sargent 1998 1 lead to the mono- 



lithic collapse model as a possible star formation scenario, 
in which every dense protostellar core collapses to a single 
star (McKee & Tan 2002 , 2003 j . However, this scenario lead 
to a time-scale problem, that can effectively destroy the sim- 



ilarity between IMF and the core mass function ( Clark et al 



20071. Similar Problems hold for the proposal that massive 



stars form by collisions of low-mass stars ( Zinnecker fc Yorke 



20071, because the observed stellar densities are too low for 
this process to be important ( Baumgardt fc Klessen|[2011 1. 
As the formation process of a cluster in a dense environ- 
ment is highly turbulent and dynamic, analytic estimates 
only give a vague idea of how accretion in this surround- 
ings may work. With the help of simulations, two different 
cases have been proposed. In one scenario, which is called 
competitive accretion, the formation of the most massive 
stars is due to a privileged position close to the centre of the 
stellar cluster, where the accretion rates are predicted to be 
highest throughout the simulation ( Bonnell et al.|[2001a|b I. 



The other scenario is fragmentation induced starvation ( Pe- 



ters et al. 2010b I, in which initially the central accretion 
rates are highest as well. However, subsequent fragmenta- 
tion shows an increasing impact of the nascent protostars 
on the accretion behaviour of the central region, that gets 
starved of material by the surrounding companions. 

In this study we investigate the accretion process in dif- 
ferent initial density profiles and different turbulent velocity 
fields. Paper I focused on the cloud structure and morphol- 
ogy. In this paper we investigate the accretion processes in 
the formed clusters in detail. 



Table 1. Physical parameters of all setups 
Parameter Value 



cloud radius 


R 


3 X 10" cm 0.097pc 


total cloud mass 


Mtot 


100 Mq 


mean mass density 


(P> 


1.76 X 10~i* gcm-3 


mean number density 


(n> 


4.60 X 10^ cm-^* 


mean molecular weight 




2.3 


temperature 


T 


20 K 


sound speed 


Cs 


0.27kms-i 


rms Mach number 


M 


3.28 - 3.64 


mean free-fall time 


tfS 


5.02 X 10* yr 


sound crossing time 


^SC 


7.10 X lO^yr 


turbulent crossing time 


itc 


1.95 - 2.16 X 10^5 yr 


Jeans length 


Aj 


9.26 X 10^ AU ft! 0.23 i?o 


Jeans volume 


V} 


1.39 X 10^1 cm? 


Jeans mass 


Mj 


1.23 Mq 



1.76 X 10 gem and the corresponding free-fall time 
gives t« = 5.02 X 10" yr. The gas with a mean molecu- 
lar weight of = 2.3 is assumed to be isothermal at a 
temperature of 20 K, yielding a constant sound speed of 
Cs — 2.68 X lO** cms~^. The Jeans length, Aj, and the corre- 
sponding Jeans mass Af, j, calculated as a sphere with diam- 
eter Aj, are Aj = 9300 AU and Mj = 1.23 Mq, respectively. 
Table [l] provides an overview of all physical parameters. 



2.2 Numerical Code 

The simulations were carried out with the astrophysical code 



FLASH (Fryxell et al. 20001 in version 2.5. To integrate 



the hydrodynamic equations, we use a piecewise-parabolic 



method (PPM) ( [Colella fc Woodward] [T984| |. The code is 
parallelised using MPI. The computational domain is sub- 
divided into blocks containing a fixed number of cells with 
an adaptive mesh refinement (AMR) technique based on the 
PARAMESH library ( [Olson et al.|1999[ ). 



2.3 Resolution and Sink Particles 

The simulations were run with a maximum effective resolu- 
tion of 4096^ grid cells, corresponding to a smallest cell size 
of Aa; ~ 13 AU. In order to avoid artificial fragmentation, 
the Jeans length has to be resolved with at least 4 grid cells 



(Truelove et al. 19971. At the highest level of refinement, 
we set the minimum resolution of the Jeans length to 6A2;, 
which can be transformed to a maximum density of 



2 NUMERICAL METHODS & INITIAL 
CONDITIONS 

The simulation data used in this study are the same as in 
Paper I, where a detailed description of the initial setups 
can be found. Here we only summarise the key parameters. 



2.1 Global Simulation Parameters 

We simulate the collapse of a spherically symmetric molec- 
ular cloud with a radius of i? = 0.1 pc and a total 
mass of 100 Mq . The resulting average density is (p) = 



pmax — 



4G(3A2;) 



= 2.46 X 10""gcm"^ 



(1) 



On all lower levels of refinement, the Jeans length is resolved 
with at least eight grid cells. If a cell exceeds this density, 
a spherical control volume with a radius of 3Aa; is inves- 
tigated for gravitational collapse indicators. If the collapse 



criteria (Federrath et al. 2010a I are fulfilled, an accreting 



Lagrangian sink particle is formed. While the total mass of 
the collapsing control volume is of the order of a Jeans mass 
at that density, the initial mass of the sink particle is only 
the mass corresponding to the overdensity above the maxi- 
mum gas density pmax. Therefore, the initial mass of the sink 
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Table 2. Numerical simulation parameters 



Parameter 




Value 


simulation box size 


^box 


0.26 pc 


smallest cell size 


Ax 


13.06 AU 


Jeans length resolution 




^ 8 (6*) cells 


max. gas density 


Pmax 


2.46 X 10~i**gcm-3 


max. number density 




6.45 X 10'' cm-3 


sink particle accretion radius 




39.17 AU 



* at highest level of refinement 

particle at the time of formation is typically significantly be- 
low the Jeans mass, and can be as low as IQ-^Mq. Table [2] 
lists the simulation and resolution parameters. 



2.4 Initial Conditions 

The following four density profiles were used: 

(i) Top-hat profile, p = const (TH) 

(ii) Rescaled Bonnor-Ebert sphere. (BE) 

(iii) Power-law profile, p cx r~^'^ (PL15) 

(iv) Power-law profile, p oc r~^" (PL20). 



A detailed description of the profiles can be found in Paper I. 

The turbulence is modelled with an initial random ve- 
locity field, originally created in Fourier space, and trans- 
formed back into real space. The power spectrum of the 
modes is given by a power-law function in wavenumber space 
(k space) with iSk oc k~^, corresponding to Burgers turbu- 
lence, consistent with the observed spectrum of interstellar 
turbulence (e.g., |Larson|[l98l| |Heyer fc Brunt||2004|. The 



velocity field is dominated by large-scale modes due to the 
steep power-law exponent, —2, with the largest mode cor- 
responding to the size of the simulation box. Concerning 
the nature of the k modes, compressive (curl-free) modes 
are distinguished from solenoidal (divergence- free) ones. The 
simulation uses three types of initial fields: purely compres- 
sive fields (c), purely solenoidal (s), and a natural (random) 
mixture (m) of both. The choice of these different turbu- 
lent fields was motivated by the strong impact of the nature 
of the modes on the cloud evolution, found by |Federrath| 
et al.| ( |2008]|2010b[ ). Notice also, that t he scale of the turbu - 
lence infiuences the result significantly ( Klessen et al.|2d00 |, 
with most realistic results corresponding to large-scale driv- 
ing adopted here (Klessen 2001 1, where we focus on decaying 



turbulence with compressive, solenoidal, and mixed modes 
are considered here. 

All setups have supersonic velocities with an rms Mach 
number A4 — «rms/cs ranging from A4 = 3.28 — 3.54 with 
an average of {A4) — 3.44. The sound crossing time through 
the entire cloud is tsc{R) = 7.10 x 10'' yr, and the time for 
gas with an average velocity of (A^) Cs to cross the cloud is 
ttc{R) = 2.06 X 10^ yr, respectively. 

We do not explicitly impose any rotation on the cloud. 
Due to the initial turbulence and randomly located velocity 
patterns in regions with different density, the total net ro- 
tational energy does not strictly vanish. However, the ratio 
Erot/\Epot\ is less than a few times 10~^, so only very little 
energy is deposited in global net rotation. 

We combine four density profiles with six different tur- 
bulent velocity fields (three different compositions of modes 



with two different random seeds each). Table [s] shows a list 
of all models. 



3 RESULTS 
3.1 Overview 

We follow the cloud collapse until 20% of the mass is ac- 
creted by sink particles. The simulation time, the number of 
formed protostars, and the mass of the most massive proto- 
star are listed in table |3] A column density plot at the end 
of each simulation is shown in Paper I (figures 4 and 5). 

The TH profile needs the longest time to form collaps- 
ing over-densities and confine 20 in sink particles. Dur- 
ing this time of about 45kyr, the turbulent motions can 
compress the gas in locally disconnected areas, leading to 
distinct subclusters of sink particles. The stronger mass con- 
centration in the centre of the BE setups and the resulting 
shorter collapse and sink particle formation time suppresses 
the formation of disconnected subclusters in favour of one 
main central cluster. The corresponding PL15 profiles show 
very similar cloud structure to the BE runs, but significantly 
different stellar properties. Due to the much stronger gas 
concentration in the centre of the cloud, all PL 15 setups 
form one sink particle very early in the simulation. This 
central sink particle accretes the surrounding gas at a high 
rate and can grow to a massive protostar before the turbu- 
lent motions eventually trigger the formation of collapsing 
filaments, which produce subsequent sink particles. Due to 
the stronger mass concentration, the clusters in the PL15 
runs are more compact. In both cases (PL15 and BE), star 
formation can proceed from inside out. The central region, 
which is much denser in the PL15 case, forms collapsing fila- 
ments earlier. As time proceeds, the turbulence can continu- 
ously compress material in outer regions, and thus form sink 
particles at larger distances from the centre. The PL20 pro- 
file only forms one single sink particle due to the very strong 
mass concentration. It forms very early in the simulation and 
accretes gas at a constant rate of about 2 x 10~^ Mq yr~^, 
extremely close to the analytical value of a highly unstable 
singular isothermal sphere (Paper I). This results in a to- 
tal simulation time of only llkyr, which is not enough for 
turbulent motions to form other filaments and further sink 
particles. 

For the analysis of the accretion behaviour in a dense 
cluster environment, we only consider the setups with mul- 
tiple sink particles. In oder to understand the accretion pro- 
cess, especially during the early phases of a cluster, one has 
to distinguish between different star formation scenarios and 
the resulting cloud morphology around the accreting proto- 
stars. Here we discuss two extreme cases, an early disc-like 
structure around a central protostellar object and an ini- 
tially filamentary structure that is not dominated by angular 
momentum. The first case is more pronounced in centrally 
concentrated density profiles where a net angular momen- 
tum with respect to the centre of mass is concentrated in a 
smaller and denser region. In contrast, initially flatter den- 
sity structures allow the formation of very massive filaments 
and subsequently collapsing regions within these filaments 
out to larger distances from the centre of mass. The angular 
momentum with respect to the centre of mass may be of 
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Table 3. List of the runs 


fiiid tlieir 


main simulation properties 








P 


turb. 


seed 


name 


-£'kin/|-£'pot| 




Sim /iff 


-^sink 






mode 








[kyr] 






[Me] 


TH 


mix 


1 


1 rl-m-i 


0.075 


48.01 


0.96 


311 


0.86 


TH 


mix 


2 


TH-m-2 


0.090 


45.46 


0.91 


429 


0.74 


BE 


com 


1 


BE-c-1 


0.058 


27.52 


0.55 


305 


0.94 


BE 


com 


2 


BE-c-2 


0.073 


27.49 


0.55 


331 


0.97 


BE 


mix 


1 


BE-m-1 


0.053 


30.05 


0.60 


195 


1.42 


BE 


mix 


2 


BE-m-2 


0.074 


31.94 


0.64 


302 


0.54 


BE 


sol 


1 


BE-s-1 


0.055 


30.93 


0.62 


234 


1.14 


BE 


sol 


2 


BE-s-2 


0.074 


35.86 


0.72 


325 


0.51 


PL15 


com 


1 


PL15-C-1 


0.056 


25.67 


0.51 


194 


8.89 


PL15 


com 


2 


PL15-C-2 


0.068 


25.82 


0.52 


161 


12.3 


PL15 


mix 


1 


PL15-m-l 


0.050 


23.77 


0.48 


1 


20.0 


PL15 


mix 


2 


PL15-m-2 


0.071 


31.10 


0.62 


308 


6.88 


PL15 


sol 


1 


PL15-S-1 


0.053 


24.85 


0.50 


1 


20.0 


PL15 


sol 


2 


PL15-S-2 


0.069 


35.96 


0.72 


422 


4.50 


PL20 


com 


1 


PL20-C-1 


0.042 


10.67 


0.21 


1 


20.0 



The table shows the setups with their turbulent modes, the random seed for the velocity field and their acronym used in further 
discussions. The ratio of kinetic to potential energy results from the applied scaling of the turbulent velocity and varies due to the 
random position of high- and low-velocity regions. The simulation time is given as total time in kyr and in units of the global free-fall 
time, tsim/%. The last two columns show the total number of sink particles A^sink and most massive particle Mmm, respectively. 



the same order as in the former case. However, due to the 
overall shallower density profile in the centre, the mass infall 
is slower and the timescale for the formation of a disc in the 
central region is larger than the timescale for turbulence to 
form filaments. 



3.2 Limiting Accretion Effects 

Before we discuss the two extreme cluster formation modes 
and their resulting accretion behaviour, we want to mention 
two effects that may limit the accretion in a cluster, namely 
the angular momentum barrier and accretion shielding by 
companion stars. In order to be accreted onto the surface 
of a protostar the gas needs to have a specific angular mo- 
mentum that is lower than the Keplerian angular momen- 
tum with respect to the centre of the potential, given by 
jKcpicr ~ V^GAIancr , where G is Newton's constant and 
Menc is the enclosed mass within radius r. In a cluster and 
especially in disc-like environments, angular momentum can 
quickly be redistributed due to local instabilities, turbulent 
motions, and strong accretion streams. We analyse the angu- 
lar momentum of the cloud with respect to the centre of the 
cluster and calculate the magnitude of the average angular 
momentum over thin spherical shells. 

For the effects of accretion shielding by surrounding 
companion stars, we estimate how much of the infalling gas 
can be accreted by secondary cluster members before reach- 
ing the centre of the cluster and becoming accreted by the 
primary cluster stars. Assuming a virialised cluster, one can 
derive the shielded fraction S using Bondi-Hoyle accretion 
( Bondi fc Hoyle|1944|[Bondi|1952[ ). The Bondi-Hoyle radius 
can be written as 



Rbh 



2GM 

c| + 



(2) 



where M is the mass of the accreting star, Ca is the speed 
of sound, and v is the relative velocity between the star and 
the immediate surrounding gas. In our analysis we inves- 



tigate the gas within a radius rsun = 100 AU around each 
protostar. An application of the Bondi-Hoyle analysis to the 
accretion in clusters and the cluster dynamics can also be 



found in Bonnell et al. (2001a I. The total shielded fraction 



of the gas falling onto a cluster with A'^ stars, an average 
stellar mass M, and a cluster radius Rci can be estimated 
by the fraction of the shielded surface area 



S = 



47riVi?; 



2 

BH 



N 



2GM 

c| -I- v'^ 



(3) 



In the case of high relative velocities between the stars and 
the gas, the Cg term can be neglected. Replacing v with the 
global virial velocity Hgiob of the cluster. 



[GMci 

"^"^ ^ y-R^ 



f GN 



(4) 



with Mc\ = NM being the cluster mass (assuming equal- 
mass stars for simplicity), yields 



S (X 



N' 



(5) 



In that case, the shielding is less efficient the more stellar 
objects there are in the cluster. Although this might sound 
counterintuitive, the assumption of virial velocities for the 
stars explains the strongly reduced Bondi-Hoyle radius with 
increasing number N . In the other extreme case of no rela- 
tive velocity of the star to the surrounding gas, the shielding 
becomes 



S ■ 



cIRci ' 



(6) 



linearly increasing with N. In our simulations we observe 
stars forming inside filaments with initially no or very small 
relative velocity between them and the surrounding gas. The 
closer to the centre the stars form or migrate during the col- 
lapse of the cloud, the more they decouple dynamically from 
the gas forming a strongly interacting A'^-body system with 
a gaseous background. We therefore expect the first case 
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Figure 1. Mass evolution of the most massive, central sink par- 
ticle in setup PL15-m-2. The upper panel shows the number of 
sink particles in the setup, the middle panel plots the accretion 
rate onto the most massive sink particle (Mmm ) and onto all sink 
particles (Mtot). The bottom panel shows the mass. 



(5* oc 1 /N) to be more important in the central region of the 
cluster, whereas the latter effect (S oc A'^) is expected to be 
dominant in the outskirts of the cluster and in the filaments. 
An analysis of the simulation is shown below, separately for 
the two extreme formation modes. 



3.3 Disc-like Accretion Mode 

In the first case we discuss the accretion scenario in a disc- 
like structure, concentrating on simulation PL15-m-2 as an 
example. A global accretion history for all protostars is 
shown in figure [l] where we plot the total number of sink 
particles as well as the accretion rate and the total mass 
confined in sink particles. The middle panel shows the total 
accretion rate onto all sink particles Mtot, the average ac- 
cretion rate per protostar (M) = Mtot(i)/Afsink(i), as well 
as the accretion rate of the most massive protostar Mmm, 
which happens to be the protostar located close to the centre 
of mass of the gas cloud and the cluster. 

During the first 7 kyr, the turbulent motions do not 
have enough time to significantly disturb the cloud. The 
anisotropies formed so far can be neglected and the assump- 
tion of a spherically symmetric collapse holds. The first sink 
particle forms in the centre of the cloud and accretes mate- 



A^sink = 1 



t = UA kyr 



iVsink = 2 



t = 14.5 kyr 



A^sink = 4 
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Figure 2. Column density plots of the central region with the 
disc around the most massive sink particle (subscript mm) in run 
PL15-m-2. The formation of secondary sink particles indicates 
the fragmentation into several objects, which quickly leads to the 
dissolution of the disc. Spiral arms develop and redirect the gas 
away from the central protostar, which gets starved of material. 
The images span roughly 4000 X 2000 AU. 
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rial in nearly free-fall. The initial density profile results in 



a linearly increasing accretion rate (see also Klessen et al, 
2000[|Sclimeja &: Klessen|2004l). At around t ~ 7kyr, the tur- 



T 1 III 



bulent motions form a central filament whose colliding arms 
concentrate angular momentum in the gas around the cen- 
tral protostar and form a disc-like object. Thus, the spher- 
ically symmetric free-fall approximation does not hold any 
longer. For the time between 7kyr < t < 14kyr the disc 
is stable and grows in mass and size by accretion from the 
filamentary arms. During that time, the angular momen- 
tum barrier prevents the accretion rate to increase further, 
resulting in a constant value of M « 5 x 10~* Mq yr~^. At 
t > 14 kyr, the disc becomes unstable, forms spiral arms and 
fragments into multiple objects. The column density plots 
in figure [2] show a time sequence of this short period. After 
the formation of other protostars, the accretion rate onto 
the central sink particle drops dramatically, while the total 
accretion rate onto all sink particles increases. The gravi- 
tational interactions between the protostars in combination 
with further infalling gas from the filament disturb the ini- 
tially disc-like structure and quickly destroy the disc. As a 
result, a more or less spherically symmetric cluster builds up. 
The average accretion rate (M) decreases by more than an 
order of magnitude after the formation of secondary stars. 
The fact that Mmm is on average several orders of magni- 
tude lower than (M) indicates that the available gas is not 
equally shared among the protostars but efficiently shielded 
from reaching the central region. This indicates that the pro- 
cess of fragmentation induced starvation , first described by 
Peters et al. (2010b I for disc-like structures, also works for 



more complex geometries. Small amounts of gas that can 
penetrate through the entire cluster and reach the centre 
can be seen as episodic accretion spikes. 

We emphasise the infiuence of different geometrical 



shapes of the cluster. Whereas in Peters et al. (2010b I the 



angular momentum vector was well defined and the star- 
vation effect was complete, the turbulent motions in our 
simulation allow for some accretion channels. The starva- 
tion effect can be understood by plotting the gas density 
around the central protostar (figure [3| and the accretion 
rate onto spherical shells around the central protostar (fig- 
ure |4| , including the net accretion onto gas shells as well 
as the accretion onto other sink particles. The gas density 
within a radius of r ~ 300 AU around the central proto- 
star first increases due to the global infall, forming a dense 
disc. Once other companions have formed, it decreases con- 
tinuously, which starves the central object. The accretion 
rate as a function of radius supports this starvation pic- 
ture. Before subsequent sink particles form, the accretion 
rate is roughly constant (M ~ 5 x 10""* yr"'^). Imme- 
diately after the disc fragments, the accretion front moves 
outwards to larger radii, resulting in significantly smaller 
values around the central protostar. As the protostars in 
the cluster as well as the gas undergo strong dynamical in- 
teractions, the net accretion rate exhibits local fluctuations. 
The curves for t 15.4 kyr are therefore averaged over sev- 
eral thousand years. Phases that include an accretion spike 
(see second panel in figure [2| lead to positive values at rgink, 
at other times the net accretion is negative, because some 
gas that is not bound to a protostar can enter and exit a 
spherical shell without being accreted onto any protostar. 
During the entire simulation time, but in particular after 
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Figure 3. Radial density profile around the central protostar in 
run PL15-m-2. The gas density first increases in the immediate 
proximity of the protostar due to the global infall. At later times, 
the surrounding companions branch off most of the gas, leading 
to a continuous decrease of the central gas density. 



the formation of subsequent protostars, the specific angular 
momentum in the disc-like object (see figure [2| as well as in 
the later formed spherically symmetric cluster is lower than 
the Keplerian value for all radii (figure [5|. Initially the ratio 
is close to unity indicating that gas motion in the disc-like 
object is affected by the angular momentum. This also leads 
to a deviation from the free-fall accretion rate for the simu- 
lation time i G (7 — 14) kyr, which would further increase in 
case of no angular momentum. After the disc becomes un- 
stable and further protostars form, the angular momentum 
is transfered efficiently resulting in a decreasing ratio. The 
low central accretion rates can therefore not be explained 
by an angular momentum barrier that inhibits the gas flow 
to smaller radii. A plot of the angular momentum of the gas 
as a function of enclosed mass for different times is shown 
in the appendix (figure [AT| . 

In addition we analyse which accretion shielding is dom- 
inant in different regions of the cloud. In figure [5] we plot 
the relative velocity of the protostars to the surrounding 
gas within a radius of rsurr = 100 AU — 2.5 racer around the 
protostar. For all times the velocities range from super-virial 
velocities in the centre of the cloud down to small relative 
velocities of the order of the speed of sound in the outskirts, 
indicating that the accretion shielding in the centre follows 
equation whereas in the outskirts, it is better described 
by equation (|6|. As expected, both cases of the accretion 
shielding effects apply. In the central region, the shielding 
therefore becomes less efficient with an increasing number 
of protostars. In contrast, the protostars in the outskirts of 
the cluster are much more efficient in shielding the accretion 
flow. In addition to the larger shielded fraction, the shield- 
ing at larger distances from the centre becomes even more 
efficient, because the protostars form along the dense fila- 
ments that channel the accretion streams. The area of the 
strongest accretion fiow then coincides with the position of 
the protostars that move with the gas flow. As a result, the 
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Figure 4. Accretion rate onto spherical shells around the central 
protostar in the cluster of setup PL15-m-2. The sink particle's ac- 
cretion radius is indicated by r^\^]^ . The plot shows several curves 
for t ^ 14.3 kyr, where there is only one sink particle (N = 1). 
Before the formation of secondary protostars, the accretion rate 
in the centre is roughly constant. Immediately after surround- 
ing companions have formed, the accretion front moves to larger 
radii and starves the central object. At later simulation times 
it ^ 15.4 kyr) the accretion rate varies, but stays very small for 
all curves. 



0.1 - 






-■ ../ = 


14.3 kyr, N = 


1-! 


« 




14.4 kyr, N = 


2 : 






14.6 kyr, N = 


4 - - 


^ 0.01 




14.7 kyr, N = 


6 




14.9 kyr, N = 


9 E 




: ./ t = 


15.0 kyr, N = 


10 : 




t = 


15.2 kyr, N = 


13 - 




t > 


15.4 kyr, > 


16 


0.001 


1 I 


1 I 1 I 


1 



0.1 



10 



[10^ AU] 



Figure 5. Ratio of specific angular momentum of the gas to the 
Keplerian value with respect to the centre of the cluster. For all 
times and all radii the ratio is smaller than unity, indicating that 
the gas is not prevented from moving inwards by angular momen- 
tum. After the formation of subsequent protostars the efficient 
angular momentum transport even lowers the ratio close to the 
centre. 



total gas mass in the centre of the cluster is very small. Fig- 
ure [7] shows the total enclosed mass (gas and protostars) as 
well as the enclosed gas mass as a function of radius. The gas 
mass is only a small fraction of the total mass and decreases 
strongly in the centre after the formation of multiple ob- 
jects. Therefore the central accretion rate is low in spite of a 
less efficient accretion shielding. The gas is already accreted 
further out before it reaches the less shielded region. 

Depending on the definition of the centre of the cluster 



100 



10 ^ 



1 r 



0.1 



~i — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — [~n 

protostars + 

^virial 



_I I 



iL 



0.01 



0.1 1 

r [10^ AU] 



10 



Figure 6. Relative velocity between the protostars and the sur- 
rounding gas within a radius of r^Mvi = 100 AU around every 
protostar as a function of distance from the centre of the cluster 
r for one snapshot at t = 25 kyr. The plots for different times look 
very similar. The relative velocity decreases for larger radii rang- 
ing from super-virial velocities in the centre down to sub-virial 
velocities of the order of the sound speed at large distances from 
the centre. 



(position of the most massive protostar, centre of mass of the 
protostars, centre of mass including protostars and gas), the 
most massive protostar is not located exactly at the centre 
of the cluster. In addition, the interaction with other pro- 
tostars leads to small displacements during the simulation. 
Neglecting the fact, that only little mass is left in the centre 
of the cluster, it could be possible that the low accretion 
rates onto the protostars in the central region arise from 
these displacements, if the protostar escapes form a coUi- 
mated accretion stream. In order to verify, that the central 
stars are really shielded, we need to analyse the motions of 
the central stars relative to the surrounding gas. We calcu- 
late the relative velocity dispersion of the gas with respect 
to the most massive protostar within a radius of 1000 AU 
around it, covering the vast majority of the entire cluster. 
We find velocity dispersions ranging from a/cs = 5.4— 10.6, 
resulting in a total time for the gas to reach the central star 
of 1.7 - 3.3 kyr after entering the sphere of 1000 AU. This 
time is only 5 — 10% of the simulated time. Given the low 
angular momentum in the centre of the cluster (figure [5|, 
the available gas has enough time to reach the most massive 
star regardless of displacements. 

3.4 Filamentary Accretion Mode 

If the formation of a cluster is not dominated by a disc-like 
structure but by filaments, the formation scenario of proto- 
stars is different (see, e.g.. 



Banerjee & Pudritz 2006). The 



initial density profiles with a fiat core need more time to form 
Jeans-unstable regions. During this time, the turbulence can 
form filaments in which protostellar condensations form next 
to each other along the densest part of the filament. The col- 
umn density plots in figure [S] show a time evolution of the 
filamentary collapse. During their formation, the protostars 
inherit the motion of the parental filament and move with 
the gas flow. With increasing proximity to other protostars, 
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Figure 7. Enclosed masses as a function of radius for the total 
mass (gas and protostars) and tlie gas mass. The gas is only 
a small fraction of the total mass in the centre and decreases 
strongly after the formation of multiple protostars. 



their attraction as an A'^-body system becomes stronger than 
the force between protostars and gas. The protostars then 
dynamically decouple from the filaments and accumulate in 
the central region in a more spherically-symmetric config- 
uration rather than in a flat or string-like structure. The 
initial filaments get dispersed in the central region, because 
the A''-body system efficiently stirs the gas. The formation of 
the first protostars is shown in the column density plots in 
figure [8] The mass accretion in this formation mode is plot- 
ted in figure |9] Note that figure [S] only covers a small time 
range at the beginning of the cluster formation, whereas fig- 
ure [9] covers the entire simulation time after the formation 
of sink particles. 

During the first ~ 6kyr after the first protostar has 
formed, the accretion rate onto the most massive protostar 
stays roughly constant. About 100 sink particles form along 
the main filaments within this time. Initially, most proto- 
stars can accrete gas from both sides of the filamentary arm, 
resulting in a high accretion rate. Despite the fact that the 
total accretion rate increases by a factor of ~ 10, the average 
accretion rate gradually drops by about one order of mag- 
nitude, because this mass fiow is distributed between ~ 100 
protostars. This is in contrast to the setup PL15-m-2, where 
the global accretion rate suddenly drops. As the setup is less 
concentrated, the Keplerian specific angular momentum of 
the gas with respect to the centre of the nascent cluster is 
significantly lower than in the PL15-m-2 case. Nonetheless, 
the ratio jgas/jKcpicr (figure [5| is even smaller than in the 
disc-like accretion mode (figure [lO| , indicating that the gas 
motion is not restricted by the angular momentum barrier. 
A plot of the angular momentum of the g cLS ELS 3, function of 
enclosed mass for different times is shown in the appendix 
(figure |A2|. 

The lower central mass concentration allows for ex- 
tended filaments and the formation of protostars at larger 
distances from the centre of the cluster. The relative ve- 
locity between the protostars and the surrounding gas (fig- 
ure |11[ ) drops significantly for larger radii, but is closer to 
the virial velocity for this setup in comparison to the disc- 



like structure. The global cluster dynamics should therefore 
correspond to the shielding relation in equation ([5|. How- 
ever, the relative velocities in the outskirts of the cluster are 
of the order of the speed of sound which marks the transi- 
tion between the two extreme shielding cases. In addition, 
the gas distribution in the filamentary structure strongly de- 
viates from spherical symmetry. The protostars in the out- 
skirts of the nascent cluster move along the densest part of 
the filament and can therefore efficiently accrete a signifi- 
cant fraction of the filament mass before they dynamically 
decouple from the filament. As a consequence, the accretion 
shielding is more efficient, the more protostars form along 
the filaments. Similar to the disc-like case, the gas content 
in comparison to the total mass is very low in the central 
region of the cluster (figure[T2|. Overall, the more protostars 
accumulate in the central region forming a cluster, the more 
efficient is the accretion shielding effect. In the case of the 
most massive sink particle in BE-m-2, accretion is entirely 
shut off for t > 27kyr. 

3.5 Mass Evolution of all Runs 

In general, the formation scenario of protostars in dense clus- 
ters will be a mixture of star formation in a disc and fila- 
ments, where the above examples are extremes. However, 
in all simulated cases, the formation of multiple protostars 
finally leads to a shielding effect in the formed cluster. Fig- 
ure [13] shows the mass evolution of the most massive proto- 
stars of all clusters combined, figures [T4||16| show the mass 
evolution of the 20 most massive sink particles for each sim- 
ulations. In the case of an initially uniform density (TH, 



figure 14 1, the two main subclusters for each run (TH-m-1 



and TH-m-2) were evaluated separately. The corresponding 
accretion plots for the most massive sink particle are very 
similar to the ones shown in figures [T] and [9] Despite different 
formation scenarios, all setups show very similar structures 
and emphasise the starvation effect on the most massive pro- 
tostars, which are located preferentially closer to the centre 
of the cluster than low-mass companions. 

Depending on the formation time, the formation loca- 
tion and the accretion rate over time, the mass of the most 
massive protostar in relation to the mass of the other ob- 
jects or the entire ensemble of protostars in the cluster might 
differ for different setups and initial conditions. Figure [17] 
shows the mass of the most massive sink particle, Mmm, in 
relation to the total mass confined in sink particles. A/tot. 
The setups with strong initial mass concentrations (PL15) 
form an early protostar that stays the only sink particle for 
a significant part of the total simulation. After the forma- 
tion of subsequent protostars, the mass of the central mas- 
sive one stays almost constant due to the starvation effect. 
The mass relations of all PL15 setups therefore bend up- 
wards. In the setups with initially flat density distribution, 
the mass of the most massive protostar roughly follows a 
relation Mmm oc M^^^ , but with remarkable scatter. The 
fact that the difference in mass between the most massive 
and the other stars is much less, leads to a stronger impact 
of A'-body interactions on the location of the most massive 
protostar within the cluster. While the very massive central 
stars in the PL15 runs remain within the accretion shielded 
area of the cluster, the most massive stars in the BE runs 
and in the subclusters of the TH setups follow larger orbits 
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Figure 14. Mass evolution of each sink particle for the four main subclusters in the TH setups. The upper two plots correspond to the 
biggest subclusters (SCI, SC2) in run TH-m-1. The lower plots (SC3, SC4) are the main subclusters in TH-m-2, respectively. The most 
massive particles are located closer to the centre of mass and thus experience an efficient starvation effect, which can be seen in the low 



increase m mass. 
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Figure 15. Same as figure [l4] but for the BE setups. 
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where they leave and reenter the shielded area alternately. 
Whenever the most massive sink particle gains further mass, 
either by leaving the accretion shielded area or experiencing 
episodic accretion, the curve flattens. Entering the accretion 
shielded volume of the cluster, the most massive object gets 
starved of material and its accretion stops. As the whole 
cluster continues to grow, the curves bend up. 



4 DISCUSSION 

Two accretion models are widely discussed in the literature: 
1) the monolithic collapse and 2) competitive accretion and 



et al. 



fragmentation induced starvation, where mutual dynamical 
interactions of accreting protostars play a key role in shaping 
the final mass distribution. The monolithic collapse model 
is based on the apparent similarity of the core mass func- 
Motte et aLl[T998l [Johnstone eral][2000l [2006] [Lada 



tion 



2002 



2008 1 and the stellar initial mass function ( Kroupa 
Chabrier||2003[ ). It rehes on a direct mapping of cores 



to stars with a constant efficiency factor for the collapse of 
a single core into a single stars or at most a binary sys- 
tem. A major caveat in this model is the collapse of very 
massive cores. The amount of turbulent energy deposited 
in these cores is likely to cause the core to fragment into 
many objects. This work and especially Paper I demonstrate 
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Figure 16. Same as figure |14| but for the PL15 setups. Note that the masses are plotted in Unear scale in order to better see the 
starvation effect, which starts after the formation of the second sink particle (vertical lines). 



that only strongly concentrated density profiles may prevent 
further fragmentation in an isothermal environment. Radi- 
ation feedback tends to reduce the degree of fragmentation, 
but still can not prevent the collapsing core from fragment- 
ing into many protostars ( Krumholz et al.|2007 Bate|2009 



[Peters et al.|[2010a[ ). Likewise magnetic fields support the 
agglomeration of larger masses. However, even magnetised 
cores fragment into smaller objects ('ZieglerpOOS Hennebelle 



fc Teyssier 2008: Banerjee et al.|20 09 Biir zle et al.|2011||Pe 



ters et al.|201lj jHennebelle et al.|2011| |Seifried et al.|2011[ ). 

The degree of fragmentation and the resulting dynamical 
accretion in our simulated dense clusters significantly differ 
from the monolithic collapse model. The fragmentation is 
suppressed only for highly concentrated density profiles and 
weak initial turbulent motions (see Paper I). Apart from 
the fragmentation of the gas, the assumption of a mono- 
lithic collapse leads to a time-scale problem that is likely 
to destroy the apparent similarity of the IMF and the core 
mass function ( [Clark et al.||2007[ [Smith et al.|2008 l 



In competitive accretion ( Bonnell et al.|[2001a|b Bon- 
nell & Bate[[2002| [Bate & BonneU[[2005[), the cloud first 



fragments down to objects with masses close to the opac- 
ity limit, which then start to accrete and build up the initial 
mass function. The accretion rates are mainly determined 
by the position of the protostar in the cluster. The gravi- 
tational potential funnels gas to the centre of the cluster, 
leading to higher accretion rates onto the most central ob- 
jects, located close to the centre of the potential well. Thus, 
central protostars can grow to the most massive ones. The 
more massive a star, the faster it can dynamically migrate to 
the centre of mass, continuously ensuring a gravitationally 
privileged position. Assuming a negligible impact of the sur- 
rounding protostars on the global mass accretion, the most 
massive stars continue to accrete at a higher accretion rate. 
As soon as the impact of further fragmentation infiuences 
the global accretion process, the central accretion rates may 
vary significantly, as it is described in the fragmentation in- 
duced starvation model. In competitive accretion, the mass 



of the most mass ive star is related to the total stellar mass 

2/3 I 

as Mnim oc Mtot ( Bonnell et al.|2004 1 . This relation was as- 
sumed to be an indicator for this accretion model ( [Krumholz 
& Bonnell 2009 1 and is in agreement with observations ( Wei- 
dner fc Kroupa[[2006[ [Weidner et al][2010[ ). However, as we 
discuss below, fragmentation induced starvation can lead to 
the same behaviour. This relation is therefore not a unique 
signpost of competitive accretion. 

Intermediate scenarios between the two extreme models 
are also reported (e.g. Peretto et al.|2006 Wang et al.|20ld l, 
in which the formation model depends on the mass of the 
star. Low-mass objects form via local collapse, while massive 
protostars and their dense proximity accumulate gas from 
the global environment. 

Recent numerical models of star cluster formation in- 
dicate that the competitive accretion model in its origi- 



nal flavour (Klessen & Burkert 2000 2001 Bonnell et al 



2001a|b| [Bonnell fc Bate[[2002[ [Bate fc Bonnell|2005[ ) needs 
reinterpretation. The suggestion of higher central accretion 
rates due to the gravitationally privileged position changes 
once the formed fragments around the centre of the collaps- 
ing area have a significant impact on the global gas accretion 
process. The infalling gas may be branched off and either be 
accreted onto these objects or dynamically redirected be- 
fore reaching the centre of the cluster. [Peters et al.[ ( [2010b[ ) 
found a strong impact of subsequent fragmentation on the 
accretion flow in a 1000 Mq rotating core including radia- 
tive feedback, but without initial turbulence. The accretion 
rate onto the central object drops significantly. Instead of 
growing further, the central protostar is starved of mate- 
rial. Despite this significantly different accretion behaviour, 
they also found masses in the simulation in good agreement 
with the relation Mmm oc M^^^ , ruling out this relation as 
evidence for competitive accretion. As they do not apply 
turbulent motions, but only solid body rotation to the ini- 
tial cloud, they focus on the starvation effect in disc-like 
conditions rather than in filamentary structures. Due to the 
missing turbulence, their disc-like structure does not evolve 
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Figure 8. Column density plots of the central filament in BE- 
m-2. The sink particles form in the filament and remain there 
while converging to the centre of the cluster. The closer the pro- 
tostars approach each other, the stronger decoupled is their mo- 
tion from the motion of the filament. The images span roughly 
4000 X 2000 AU. 
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Figure 9. Same as figure [l] but for the simulation BE-m-2 (fila- 
mentary accretion). 
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Figure 10. Same as figure |5] but for setup BE-m-2 (filamentary 
accretion). 



into a volume-filling cluster, but remains a flattened struc- 
ture. A disc provides a smaller effective area through whicii 
gas can flow towards the centre. In a cluster with random 
protostellar trajectories, the effective accretion area is much 
larger. Therefore, from simple geometrical arguments, it is 



much more difficult to shield the central area from accretion 
flows. 

However, we find that protostars in the outer region 
of the cluster form along the densest part of extended fila- 
ments, which allows them to accrete much more mass than 
in an idealised spherically symmetric setup. In addition to 
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Figure 11. Same as figure [6] but for setup BE-m-2 (filamentary 
accretion) at t = 28 kyr. 
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Figure 12. Same as figure[7]but for BE-m-2 (filamentary accre- 
tion). 



that, we observe significantly smaller relative velocities be- 
tween the protostars and the immediate surrounding gas for 
protostars at larger distances from the centre of the clus- 
ter. This increases their geometrical accretion cross section, 
expressed by the Bondi-Hoyle accretion radius. Altogether, 
the accretion in the central region of the cluster is influenced 
by efficiently shielding protostars located at the position of 
the accretion flow. As a result, even a relatively small num- 
ber of surrounding stars is sufBcient to provide an efficient 
starvation effect in a volume-filling cluster environment. 

One caveat of the presented simulations is the missing 
feedback from the stars that are formed first in the clusters. 
The gas stays isothermal and therefore tends to fragment 
much more quickly than with heating of the gas. Includ- 
ing this effect, the number of protostars is expected to be 
smaller. The starvation effect on the central objects in the 
protostellar cluster might thus be overestimated, especially 
in a turbulent volume-filling cluster, where gas can fall into 
the central region from all directions. However, [Peters et al.| 
1 2010b I found that the number of protostars in simulations 
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Figure 13. Mass evolution of the most massive protostar for all 
simulated clusters. In all cases the mass curves flatten significantly 
at the late-time evolution of the simulation. For a detailed mass 
analysis of each setup, see figures [T4j]16| 
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Figure 17. Relation of the most massive star in the cluster to 
the total cluster mass. Note that in the PL15 case, secondary 
sink particles form very late in the simulation (indicated by the 
arrows), which explains the large range in mass with a slope of 
unity. 



including radiative feedback is roughly half of the number of 
protostars in an isothermal calculation. As the total number 
of sink particles in our simulations is much larger than 100, 
and additionally, the starvation effect can already be seen 
with a relatively small number of competing protostars, we 
do not expect the shielding effect to completely vanish in 
the non-isothermal case. As long as there is some fragmen- 
tation, the fragments are likely to shield the central massive 
star from accretion, irrespective of whether radiation feed- 
back is included or not. Likewise do magnetic fields tend to 
reduce the degree of fragmentation, but still do not prevent 
the cloud from fragmenting ( Hennebelle et ar]|2011 Peters 
et aL|2011| ) 



Finally, we want to point to recent studies by |Kruijssen| 
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et al. ( 2011 1. They analysed the substructure within clusters 



as well as the dynamical state of the stellar cluster when gas 
expulsion becomes important, i.e., at a slightly later stage 
of the evolution of the cluster. Analysing the simulations of 



Bonnell et al. ( 2003 2008 1, they find that the stellar system 



quickly reaches a globally virialised state if the gas poten- 
tial is excluded and the stellar system is followed with pure 
A'^-body dynamics. Their results support the evolutionary 
picture of the formation of protostars that we see in our 
simulations. New protostars that form at larger radii from 
the centre of the cluster in gas-dominated regions have sub- 
virial velocities. As soon as they decouple from the gas mo- 
tion and move to the central gas-poor environment, they 
quickly virialise. 



5 SUMMARY AND CONCLUSION 

We performed simulations of collapsing molecular gas clouds 
with a total mass of 100 and a diameter of 0.2 pc. We 
varied the initial density profile as well as the turbulent su- 
personic velocity field and analysed the fragmentation pro- 
cess and the accretion onto nascent protostars. Most clouds, 
especially the ones with initially less centrally concentrated 
density profiles undergo fragmentation and form compact 
stellar clusters. We studied the accretion process in these 
clusters with the following conclusions: 

• All clusters show strong and fast dynamical interac- 
tions between the protostars. During the formation time of 
the cluster, the protostars complete several orbits in the 
parental cluster. The mass accretion process can thus not 
be described by a monolithic collapse model. 

• Fragmentation and formation of multiple protostars 
strongly infiuence the subsequent accretion flow in the entire 
cluster. Gas that falls towards the centre of the cluster is ef- 
ficiently accreted onto protostars that are located at larger 
radii from the cluster centre. As a consequence, the cen- 
tral region is effectively shielded from further accretion and 
none of the central objects can sustain its initially high ac- 
cretion rate. In all cases the observed decrease of the central 
accretion rate is due to the efficient shielding by secondary 
protostars. A significant starvation effect due to angular mo- 
mentum and a resulting radial barrier can be excluded. 

• In setups with initial fiat density distributions, the mass 
of the most massive protostar, Mmm, scales with the total 
cluster mass, Mtot, like M^im oc M^^^ , as originally pro- 
posed as a probe for competitive accretion. This relation is 
a common feature to both dynamical accretion models and 
cannot be used to distinguish between competitive accretion 
and fragmentation induced starvation, a conclusion that was 



already reached by Peters et al. (2010c I 



• The accretion process does not differ between density 
profiles that form only one main cluster (power-law profiles, 
Bonnor-Ebert density distribution) and setups that form 
multiple subclusters (uniform density profile). In both, the 
central clusters as well as the disconnected subclusters, we 
find fragmentation induced starvation to work in the central 
region of the cluster. 
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Figure Al. Specific angular momentum of tiie g function 
of enclosed mass for PL15-m-2. The plot shows temporal changes 
of the angular momentum of the gas during the formation of the 
cluster. 
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Figure A2. Same as figure [Al] but for setup BE-m-2. 



APPENDIX A: ANGULAR MOMENTUM 

In order to see the time evolution of the angular momentum, 
we show the specific angular momentum as a function of en- 
closed mass for the setups PL15-m-2 and BE-m-2 (figures[AT 



and A2|. In both cases the angular momentum distribution 
is exposed to variations of the order of a few up to almost 
an order of magnitude due to the turbulent interactions and 
the gas accretion onto protostars. Due to the unstructured 
motions of the protostars and the random character of the 
turbulence, the variations do not show systematic changes 
over time, i.e., the gas in the cluster can quickly gain angu- 
lar momentum due to accretion streams or lose it via shocks 
and dissipation. 



